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ABSTRACT 

Starting from the problem of two fixed centres we find a simple derivation of its third 
integral in terms of the scalar product of the angular momenta about the two fixed 
centres. This is then generalised to find the general form of the potential in which an 
exact third integral exists. 
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1 INTRODUCTION 

When deriving or using the third integral of the motion, 
i.e., that other than energy or angular momentum about 
the axis, I have often been asked for its physical meaning. 
The new derivation given here may not completely answer 
that question, but the realisation that the kinetic part is the 
dot product of the angular momenta about the two centres, 
or two foci, goes somewhat further towards an answer than 
anything I have seen previously. 

This paper gives what I hope is a simpler, cleaner 
derivation of the integral, in elementary terms that do 
not require spheroidal coordinates or the separation of 
the Hamilton-Jacobi equation. Earlie r treatments requir - 
ing such s o phisti c ations are found in Lyndcn-Bcll ( 1962 ), 
Eddingtoij (jl91S| ), |de Zeeuw| ( |l985abd) 7 ptackcl| fll89C| ) and 
Kuzmin _[ 1956 ) . T he two dimens i onal p roblem is discussed in 
Whittakei (1904). |Contopoulos (1960) showed that in many 
galactic potentials of non separable form, a regular third 
integral exists for most orbits. 



2 PRELIMINARIES - THE ONE CENTRE 
PROBLEM 

Hamilton's beau tiful treatment o f this problem is not widely 
enough taught (Hamilton, 1847). As we need some of the 



steps later we give his treatment here. The equation of mo- 
tion is 



-fir/r 



(1) 



where fi — Gm and f is the unit vector. Cross multiplying 
by r and writing h = r x r we find 

h = r x r = , 

so h is constant. Cross multiplying (lit) by h 



d(hxf) jdt — —fx (f x r/r) x f = —fxdr/dt 



(2) 



the last identity may readily be demonstrated by dull alge- 
bra but the physical understanding is more interesting. In a 
small interval St, vx8v/r is perpendicular to the movement 
<5r and to f and gives the angle moved about the centre. 
Hence f x r/r is the angular velocity of the radius vector 
which we call uj. Hence (f X r/r) X f = u X f where uj is the 
angular velocity of the unit vector f . However, the change 
in a unit vector is just uj x f because it only changes by 
rotation since its length is fixed so dr/dt — uj x f . 
Equation (JiJ) is readily integrated to give 



r x h = fi(r + e) 



(3) 



where e is a constant of integration. Since both f and f x h 
are perpendicular to h we find e ■ h = . Take the scalar 
product of (^) with r/fi then, on use of 



r • (r x h) = (r x r) • h/r — h /r 
we find 



£/r = l + e- 



(4) 



where I — h 2 /fi. This is the equation of a general conic 
where e ■ f = e cos <j> with <j> measured from pericentre and e 
the eccentricity. So we now identify Hamilton's vector e as 
the vector that points to the pericentre and whose length is 
the eccentricity. We note in passing the useful expression for 
the velocity found from (H), 



r = /ih h x (f + e) . 

Its square may be used to derive the energy integral 



(•>) 
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3 THE THIRD INTEGRAL FOR THE TWO 
CENTRE PROBLEM 

We measure ri from the first centre of attraction at —a and 
r 2 from the second at +a so ri = r 2 + 2a. The equation of 
motion is 



^d/dr 2 (rx^tp) + 



ri = r 2 = -nih/rj - liii^lri 



(6) 



The energy and the component of angular momentum along 
the line of centres are, of course, conserved. We are interested 
in an independent third integral. Cross multiplying (^) by 
ri we find 

hi = d(n x ri) jdt = -/wri x r 2 /V 2 = -^ 2 2a x r 2 /r 2 , 

taking the scalar product with h 2 we find c.f. (^), 

h 2 ■ hi = — fj, 2 2a ■ [i 2 x (f 2 x r 2 /r 2 )] — /i 2 2a ■ dr 2 /dt . 

Adding a half of this to the conjugate expression with r 2 
and ri exchanged and —a for a we find 

d ( i hi ■ h 2 ) fdt = d[a ■ (fi 2 r 2 - fiih)] fdt . (7) 

So 

h = jhi ■ h 2 + (/wri — ji 2 i 2 ) ■ a = constant , (8) 
which is the third integral for this problem. Using 
r=|(n+r 2 ) ; ajl = (r± - r 2 )/2 

as coordinates the three integrals may now be used to solve 
for the orbits, but our purpose here is to generalise this 
beautifully simple derivation of to more general prob- 
lems. When \x 2 = 0, h 2 — hi — 2a x ri and the integral 
reduces to a combination of hi and e. 



4 THE EXACT THIRD INTEGRAL MORE 
GENERALLY 

We write our axially symmetrical potential ip in the form 
■0(n, r 2 ). The equation of motion is 

f i = r 2 = Vy> = dtp/driii + dip/dr 2 r 2 . 

Performing the same steps as led to (Q) we find in its place 

d(ihi ■ h 2 ) jdt = —r 2 dip/dr 2 a-dr 2 /dt+r 2 dip/dria-dri/dt , (9) 

which reduces to our former result when ij) = y,\/r\ + fi 2 /r 2 . 

To proceed we need expressions for a ■ r 2 , etc. These 
come from squaring the expression for ri in terms of a 
and r 2 . 

r'i=r 2 + 4a 2 + 4a ■ f 2 r 2 



r 2 d(a- f 2 ) jdt = \r\r 2 r\ 



(n + r 2 ) 



a \r 2 



(10) 



Collecting terms in ri and in r 2 in (|9() we have using ( ]1 0| 
and its conjugate with n < > r 2 and a — > —a 



d (ihi ■ h 2 ) jdt = {-lr ir2 di>/dr 2 ) 
+ [i (n + r 2 ) - a 2 ] dtp/dri} fi + conjugate . 
This last expression may be rewritten 



(11) 



+ id/dn [(r 2 + r|-4a 2 )v] , 



> ri + conjugate 



The second term in the curly bracket gives rise to an expres- 
sion which, taken with its conjugate is a total time deriva- 
tive. 

We now ask for the condition on ip that the remainder 
plus its conjugate be a perfect derivative. For this to be true 
it must be of the form dA/dri fj + dA/dr 2 r 2 so d/dr 2 of 
the first term must be d/dri of its conjugate. 

Writing \ = rir 2 ip our condition reduces to 

d 2 x/drl - d 2 x/drl = , 

but this is the wave equation and its general solution is 
X = C(Ti + r 2 ) - r)(n - r 2 ) , 

where ( and -q are arbitrary functions of those arguments. 
Thus the potential tp must take the form 

ip = [((ri + r 2 ) - r)(n - r 2 )] /{r\r 2 ) . 

With this expression for ( p^ ) takes the form 

d 



dt 



(ihx-ha) 



dri 



+conjugate, 



so 

h 



3 hi ■ h 2 
ihi ■ h 2 
1 

rir 2 



[i {r\ + r 2 ) - a 



r\r 2 



+ hiC + V) = const 



n — r 2 



c- 



n + r 2 



The final term reduces with suitable definitions to the 
— (nC — XrfjJjX — /i ) form one finds in standard treatments 
( |Lynden-Beli| |l962j ). 

When the coordinates are oblate a is imaginary i|a| so 
hi and h 2 are angular momenta (r + i|a|) x r, etc, about 
the imaginary points r — ±i|a|. However the product of 
these angular momenta is still real. The expressions for the 
integral are unchanged except that when the r\ and r 2 are 
interpreted as the greatest and least distances of a general 
point from the ring z = 0, r = \a\ then the —a 2 terms in 
the final expression of 73 should be omitted. 

I am unable to extend this derivation into the relativis- 



tic regime while maintaining its simplicity. Carter ( 196S ) 
has shown that a third integral exists for a charged particle 
moving in the Kerr-Newman metric. When Newton's G is 
set equal to zero that metric leaves behind the most inter- 
es ting electrom a gneti c field with a net charge qi (discussed 
in [Lynden-Belil (|200l[)). 



E-HB = -V qi/yj{r-ia) ■ (r - ia) 



(0,0, a) 



Third integrals in this and other electromagnetic fields in 
the relativistic regime are related to the classical third inte- 
gral discussed above. For a demonstration of this relation- 
ship and a derivation of the forms o f elect romagnetic fields 
in which they exist see Lynden-Bel] ( 2000 ) . Here we merely 
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quote the form the third integral takes for the above elec- 
tromagnetic field with the orbiting particle of mass m and 
charge q. is the unit toroidal vector, $ and A = Acf> the 
electromagnetic scalar and vector potentials. 



I = 


Jhi •h 2 + gK3(A)/(A-M) 


where 




hi = 


(r — ia) x p ; I12 = (r + ia) x p 


Cs(A) = 


qi \/ X — a 2 A _1 c~ 2 ^eA — ^qqi \J X — a 2 — ach 


P = 


mv/y/l — V 2 /c 2 + qAcj)/c ; v = r 


e = 


mc 2 /^l-V 2 /c 2 + q<S> 


h = 


mR 2 <j>/ y/l - V 2 /c 2 + qRA/c 


$ = 


qi\/X- a 2 /(X ~ fi) 


A/R = 


a 2 A _1 $ 


2A = 


r 2 + a 2 + (r — ia) 2 



(constant on oblate confocal spheroids) 

2/1 = r 2 + a 2 — I (r — ia) 2 | 

(constant on confocal hyperboloids of one 
sheet) 

7, h and e are the integrals, R 2 — x 2 + y 2 = Xfx/a 2 , and the 
I sign is used in the sense of complex numbers. 
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